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We present further results on a class of sums which involve complex powers of the 
distance to points in a two-dimensional square lattice and trigonometric functions 
of their angle, supplementing those in a previous paper (McPhedran et al, 2008). 
We give a general expression which permits numerical evaluation of members of the 
class of sums to arbitrary order. We use this to illustrate numerically the properties 
of trajectories along which the real and imaginary parts of the sums are zero, and 
we show results for the first two of a particular set of angular sums which indicate 
their density of zeros on the critical line of the complex exponent is the same as 
that for the product of the Riemann zeta function and the Catalan beta function. 
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1. Introduction 

This paper adds to results in McPhedran et al (2008) (hereafter referred to as I) 
on the properties of a class of sums over two-dimensional lattices involving trigono- 
metric functions of the angle to points in the lattice, and a complex power 2s of 
their distance from the lattice origin. There, it was shown that certain of these 
angular sums had zeros on the critical line Re(s) = 1/2, but could not have zeros 
in a neighbourhood of it. 

We derive a general expression which is exponentially convergent and permits 
the rapid and accurate evaluation of the angular sums irrespective of the value 
of the complex parameter s. We demonstrate the high-order convergence of this 
formula by using it to illustrate a limiting formula for a particular set of angular 
lattice sums. We go on to consider the properties of trajectories along which the 
real and imaginary parts of a class of angular sums are zero, and in particular we 
establish accurate approximations for these trajectories when Re(s) lies well outside 
the critical strip < Rc(s) < 1. We give preliminary results on the distribution of 
zeros on the critical line Re(s) = 1/2 of two angular sums, which suggest that to 
leading order they have the same density of zeros as the product of the Riemann 
zeta function and the Catalan beta function. In Appendices we comment on the 
functional equation satisfied by a class of angular sums, and the properties which 
link angular sums of order up to ten. 

The analytic results presented here are supported by numerical results obtained 
using Mathematica 7.0.1. Formal proofs of certain key properties of the angular 
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sums and the location of their zeros will be given here, while others will be presented 
in a companion paper. 

There are two principal motivations for the study presented here. The first 
is that the very general expression for the class of angular lattice sums derived in 
Section 2, and their connections with other angular lattice sums shown in Appendix 
A, enables them to be used in physical applications requiring regularization of 
sums over a two dimensional square lattice. The class of summands which can 
be addressed is wide, as it consists of any function which has a Taylor series in 
integer powers of trigonometric functions of the lattice point angle in the plane, 
and any complex power of its distance from the origin. For example, any summand 
of the type often encountered in solid state physics combining a Bloch-type phase 
factor and a function of distance having a Taylor series could be so represented. 
The second is that we show the angular lattice sums to be connected with the 
product of the Riemann zeta function and the Catalan beta function in a very 
natural way- for example, it seems that the densities of their zeros on the critical 
line are the same to leading orders. In this way, these angular sums may provide 
a new way to forge a link between the Riemann hypothesis and its generalization 
to other Dirichlet L functions, as well as providing a wide class of functions, in 
which identified members obey a priori the hypothesis, and others do not. There 
are interesting parallels between this work and that of S. Gonck (2007), although 
we deal with double sums and Gonek with single sums. 

2. An absolutely-convergent expression for angular lattice 



We recall the definition from (I) of two sets of angular lattice sums for the square 
array: 



where 9 Pl , P2 = arg(pi + i?>2), and the prime denotes the exclusion of the point at 
the origin. The sum independent of the angle PllP2 was evaluated by Lorenz (1871) 
and Hardy(1920) in terms of the product of Dirichlet L functions: 



C(0, m; s) = 5(0, m; s) = C(0, 1; s) = 4Zi(s)Z_ 4 (s) = 4C(s)L_ 4 (s). (2.2) 



Here £i(s) is more commonly referred to as the Riemann zeta function, and L-^(s) 
as the Catalan beta function. A useful account of the properties of Dirichlet L 
functions has been given by Zucker & Robertson (1976). 

It is convenient to use a subset of the angular sums (|2.ip as a basis for numerical 
evaluations. We note that the sums C(n, 1; s) are zero if n is odd. We next derive 
the following relationship for the non-zero sums C(2n, 1; s): 



sums 




(2.1) 





(2.3) 
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where K v (z) denotes the modified Bessel function of the second kind, or Macdonald 
function, with order v and argument z. The general form (|2-3[) may be derived 
following Kober (1936) in the usual way: a Mellin transform is used to give 

2n 2n poo 

Y ' , 2 2, T = y / f+n-h-^+^dt. (2.4) 

,^ , (Pl+P 2 2 ) s+n ( ^ T(s + n)J a V ' 

(Pl,P2) (Pl,P2) 

The Poisson summation formula is then used to transform the sum over p 2 , giving 

Pi™ _ / Pi™ f°° ,s+ri-l„-tp? F^„-tt 2 pIi 



Y ' vi = y , Pi t s+n-i M e -^/t dt . (2.5) 

(Pl,P2) (Pl,P2) 

We then separate the axial contribution, which for n =/= comes from p 2 = alone, 
and use Hobson's integral 

t^e-^-^dt = 2 Q 7 X S (2V?P) (2.6) 



on the remaining double sum. This leads directly to 

It should be noted that the double sum in (|2.3[) is exponentially convergent. 
Indeed, from relation 9.7.2 in Abramowitz and Stegun (1972), the large argument 
approximation for the Macdonald function of order v is 

This means that the double sum starts to converge rapidly as soon as the argument 
2irpip2 everywhere exceeds the modulus of the order s + n— 1/2. In practice, accu- 
rate answers are achieved when sums are carried out over p\ and pi from 1 to P, 
where P ~ \s + n — 1/2\/tt (the precise value of P required being fixed by studies 
of the effect of increasing P on the stability of the result). The representation (|2.3[) 
and finite combinations of it thus furnish absolutely convergent representations of 
trigonometric sums from the family (|2.1[) and close relatives, for any values of s 
with finite modulus. These representations are easily represented numerically in 
any computational system incorporating routines for the Riemann zeta function of 
complex argument, and Macdonald functions of complex order and real argument. 
(In practice, the Macdonald function evaluations are most time-expensive in the re- 
gion of (pi,p 2 ) values where pip 2 has comparable magnitude to \s + n — l/2\/(2n). 
Thus, it is efficient to create a table of these values for p\ varying with p 2 — 1, run- 
ning up to an argument where the Macdonald function is less than an appropriate 
tolerance times its value for p\ = 1.) 

As an example of the numerical efficacy of (|2.ip , we consider its use in illustrating 
a limiting property of the sums C(2m, 1; s).We have 

(Pl,J>2) 

and as m — > oo we require | cos9 pi P2 \ — 1 for a contribution, i.e. 

lim C(2m,l;s) = 2((2s). (2.9) 
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Figure 1. The modulus of C(2m, 1; s) (red) and 2^ (2s) (green) as a function of 
s = 1/2 + it, for t £ [0, 20], with m = 10 (left) and m = 100 (right). 

The relationship (|2.9[) is illustrated numerically in Fig. Q] For the right-hand side of 
(|2.9p to be accurate, the required order m increases with t, although convergence 
is also slow for t near zero. 

Another angular sum of great importance in this paper can easily be expanded 
in terms of the C(2n, 1; s): 



C(1,4to;s)= J2 



(2.10) 



or, in terms of the Chebyshev polynomial of the first kind (Abramowitz & Stegun 
(1972), Chapter 22), 



C(l,4m;«) = 



Pl,P2 



i T^m (cos 9p 1 t p 2 ) 



(2.11) 



As the coefficients of this Chebyshev polynomial are explicitly known, the represen- 
tation (|2.11[) enables any sum C(l,4m;s) to be expressed as a linear combination 
of sums C(2n, 1; s) with < n < 2m. 

The connections between various angular lattice sums grouped in systems with 
order up to 10 are explored in Appendix A. 



3. Some properties of trigonometric lattice sums 

The functional equation is known (see McPhedran et al. (2004), eqs. 32 and 59) for 
C(l,4m;s): 

Pi's J- Om \ 

G im {s) = C(l, 4m; s) [S+ g ' = G im (\ - s). (3.1) 

This equation also holds for m = 0, where it gives the functional equation for 
the product C(s)L_4(s). It is in fact the m dependence of the functional equation 
(|3.1jl which enables the derivation of many of the results in (I) and the present 
paper. As the derivation in McPhedran et al (2004) uses different notation to that 
in subsequent papers and here, we give a brief discussion of the argument leading 
to (|3.1[) in Appendix B. 

This m dependence in (|3.ip is represented in two related functions: 

^- / r(l — s + 2m)T(s) . , , . . , „. 
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where 4>2m{s) is in general complex. Note that J-2 m {s) is the ratio of two polynomials 
of degree 2m, with one obtained from the other by replacing s by 1 — s: 

(2m- 8 )(2m-l-s)...(l-s) 
2m{ ' ~ [2m - (1 - s)][2m - 1 - (1 - a)] . . . [1 - (1 - a)] ' [ ' 

We then introduce two sets of rescaled lattice sums: 

£7(2, 2m; a) = — -^=L={C(0, 1; s) + £7(1, 4m; a)}, (3.4) 

27T s V-^2m(s) 

and 

5(2, 2m; a) = - ^== [C(0, 1; a) - £7(1, 4m; *)]. (3.5) 

27T S V /J r 2m(s) 

Note that this definition means C and S 1 have branch cuts where ^2m(s) is real and 
negative. For example, for .7-2 (s), the branch cut includes a circle in the (er, t) plane 
(s = g + it), with centre (1/2,0) and radius V3/2. 

An analytic function which combines lattice sums is defined by 

A 3 (2, 2m; a) = C(2, 2m; s) 2 - 5(2, 2m; s) 2 = C(0, 1; s)C(l, 4m; s). (3.6) 

It obeys the functional equation 

A 3 (2, 2m; s) = T 2rn {\ ~ s)A 3 (2, 2m; 1 - s), (3.7) 
and on the critical line s = 1/2 + it, 

A 3 (2, 2m; a) = [1 - itan(<fe m (s))][|C(2, 2m; s)| 2 - |5(2, 2m; s)\ 2 ], (3.8) 

and 

Im[C(2, 2m; 1/2 + it)} = Im[5(2, 2m; 1/2 + it)]. (3.9) 
We note from $3M that 



ImA 3 (2,2m;i+it) = - tan(0 2m , c t) Re A 3 (2, 2m; i + it), (3.10) 

using the notation (1/2 + it) = </>2m,c(t), a real- valued function. We can take 
the derivative with respect to t of (|3.10p , to obtain: 

-^ImA 3 (2,2m;~ + it) + tan 2m , c (t)^ Re A 3 (2, 2m; ^ + it) 

COS^(02m,c(i)) 2 

We use the Cauchy-Riemann equations in (|3.11[) . to obtain 

^ReA 3 (2,2m;i + it) + tan 2m , c (t)^ Re A 3 (2, 2m; ^ + it) 

ReA 3 (2,2m;i+it), (3.12) 



J (02m,c(t)) 
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and 



5 1 d 1 

— ImA 3 (2,2ro;- + it) - tan0 2m , c (t)— Im A 3 (2, 2m; - + it) 

at 2 oo~ 2 

02m,c(*) . ReA3 ( 2j 2TO;i+ii). (3.13) 



COS 2 (0 2m , c (t)) " '2 

The equation Q3.12[) indicates that, at points where contours of Re A3 (2, 2m; a + 
it) = intersect the critical line, their tangent vector is given by (1, tan</>2m,c(^))) 
provided 9 A3 (2, 2m; | +it)/dt 7^ 0. (This requirement is that the left-hand side 
of (|3.12[) can be interpreted as the scalar product of the tangent vector and the 
gradient of Re A3, with the latter having a well-defined direction.) The equation 
(|3.13p indicates that the tangent vectors for the contours of Im A 3 = at the critical 
line are given by (— tan 02m, c(t), l),i.e. they are at right angles to those for the real 
part. 

We note that, for \s\ » 1, or, more strictly, for |s| >> 4m 2 , 

4m 2 8m 4 2m 2 i im 2 (8m 2 — 1) 

T 2m { s )^l-— Yf2+{s _ i/2r hmis )^-— J -+ 6(s _ 1/2)3 . (3.14) 

Thus, for \t\ >> 1, (|3.12p takes the approximate form 

d 1 2m 2 d 1 2m 2 1 

— Re A 3 (2, 2m; - + it) + — — Re A 3 (2, 2m; - + it) = _ R e A 3 (2, 2m; - + it). 
Oct 2 t at 2 t z 2 

(3.15) 

This shows that, as t increases, the contours of Re A^(2, 2m; a + it) = strike the 
critical line at ever flatter angles (although the angle increases as m increases) . The 
direction of the intersection with the critical line is unique for every point where 
<9A 3 (2, 2m; i + it)/dt ^ 0. Similar remarks apply to Im A 3 (2, 2m; er + it) = 0, where 
the equivalue contours cut the critical line at a direction given by the tangent vector 
(— 2m 2 /t, 1), and thus their gradient at the point of intersection increases with t. 



4. Equivalue contours of Re A 3 and Im A 3 

In Figs. [2E1 we show contours on which the real and imaginary parts of A3 (2, 2; s) 
are zero. Fig. [2] gives some detail of the region near a — 0.5 for t e [0.1, 10], and a 
further more detailed region on an undistorted geometric scale (including the semi- 
circle on which Ti[s) is real). Fig. Ogives a more global view of null contours for t 
ranging up to 20. Fig. [4] gives the detail of the null contours for t near two values at 
which contours nearly touch. We introduce the notation 14 for a zero on a — 1/2 of 
C(l, 4; s), +1 for a zero of C(s) and —4 for a zero of L_4(s). Then the zeros evident 
in Fig[3] are categorized as: -4, 14, 14, -4, 14, -4, +1, 14, -4, 14, -4, 14, 14. 

The null contours have a number of interesting properties. Firstly, we can see 
that contours for the real part do indeed strike the critical line at a small angle 
to the a axis, which decreases as t increases, while the contours for the imaginary 
part strike the critical line almost vertically. Secondly, the contours for the real and 
imaginary parts intersect the critical line simultaneously, except for one point. This 
is at (1/2, v / 3/2), where, as remarked in (I), 02(1/2 + 1^/3/2) = n/2, which permits 
ReA3(2,2;s) to be zero, while ImA 3 (2,2;s) is non-zero (see (|3.8|) V 
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Figure 2. Null contours of the real part (red) and imaginary part (blue) of A3 (2, 2; a + it), 
with (left)cr £ [—0.5, 1.5], and t £ [0, 10]. On the right we show part of the left panel, with 
now t £ [0.0,2.0], to give an undistorted geometric image. On the green circle, J-2(s) is 
real and negative. 

We can see four null contours for the real part and five null contours for the 
imaginary part intersecting the axis t — 0. If we use the expression 



A 3 (2,2; S ) 



r(s)a(a + l) 

(2-s)(l-s)7T 5 



-[C 2 (2,2; S )-S 2 (2,2; S )], 



(4.1) 



we find that A 3 (2,2;s) has a second order pole at s = 1, and first order poles at 
s = and s = 2. Numerical investigations show that near these points 



A 3 (2,2;1 + <S) 



4.59643 



and 



A 3 (2,2;<5) 



-0.798212 



A 3 (2,2;2 + <5) 



1.16981 



(4.2) 



(4.3) 



6 ^ov,,.-, g 

We note that ReA 3 (2,2;s) has zeros at the two first-order poles, where the tra- 
jectories approach the poles broadside (i.e., parallel to the t axis), and two null 
trajectories approaching the second-order pole at 45° and 135° to the a axis. 
ImA 3 (2,2;s) has a zero at the second-order pole, where the null trajectory ap- 
proaches broadside. The other null trajectory crossings of the real axis occur at 
0.29782, 1.67735, -2.65568 and 4.21422. The last two of these are associated with 
minima of Re A 3 (2, 2; a). 

We show in Fig. the behaviour of the null contours for a large enough (a > 4) 
to enable us to accurately approximate trigonometric sums by their first few terms: 



C(2,2; S ) 2 -5(2,2; S ) 



16(1 



1 



36 

52S+: 



so that 



A 3 (2,2;s) 



i6r(i-s)r(s + 2)r(s) 



r(3-s)7r 



2 s 



(1 + — 

V 4 s 



), 



36 , 



(4.4) 



(4.5) 
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5 10 



Figure 3. Null contours of the real part (red) and imaginary part (blue) of 
A 3 (2, 2; a + it), with a £ [-4.5, 11.5], and t G [0.1, 20]. 

For a beyond 4.21422, we see from Fig. [5] that neither the null trajectories of the 
real part or the imaginary part of A 3 (2,2;s) can attain the real axis. In Fig. [SJ 
we show for comparison the null trajectories associated with the prefactor term in 
square brackets in (|4.5[) . It is evident from a comparison of Figs. [5] and [5] that the 
null trajectories are given accurately by the prefactor in the range of a shown. Using 
Stirling's formula (Abramowitz and Stegun (1972)), we can place the prefactor in 
exponential form for |s| large, and find the constraint for null trajectories to be 

s 25 2 1 

Im[2s log(-) - 2s - log(s) + — + —] = (n + -)tt, or nn, (4.6) 

7T OS S 2 

where on the right-hand side the first value is for real-part zeros, and the second 
for imaginary-part zeros (n being an integer). From (|4.6[) we find that the equiv- 
alue contours alternate for a large, with null contours for the real part sandwiched 
between those for the imaginary part, and vice versa (with each trajectory corre- 
sponding to a first-order zero). Each contour tends to zero as 1/ log(cr/7r) as a — > oo. 
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-10 12 -10 12 



Figure 4. Detail of the null contours of the real part (red) and imaginary part (blue) of 
A 3 (2, 2; a + it), for t near 10 (left) and near 18.5 (right). 




4 5 6 7 

Figure 5. Equivalue contours of the real part (red) and imaginary part (blue) of the form 
g2) for A 3 (2, 2; a + it), with a G [3.5, 7.5], and t € [0.1, 10]. 

Note that we may use the equation (13. 7|) to deduce relations governing the phase 
of A3(2, 2m; s). Dividing (|3.7|) by its conjugated form, we obtain 

arg(A3(2,2m; S )) + arg(A 3 (2,2m;l-s)) = - arg(JF 2m (s)). (4.7) 

This agrees with (13.81) and (|3.10[) when Re(s) = 1/2. Furthermore, on the right-hand 
side of Fig. [3] we know that arg(A 3 (2, 2; s) increases monotonically with t, from its 
value of zero on the real axis. This enables us to assign values of the constant phase 
of the null contours there: tt/2, it, —it/2, 0, ir/2, etc. Using (|4.7p and (|3. 14[) . we see 
that the monotonic increase of arg(Aa(2, 2; s) with t on the right in Fig. [3|) forces 
a monotonic decrease of arg(Aa(2, 2; s) with t on the left. We can then assign the 
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4 5 6 7 



Figure 6. Null contours of the real part (red) and imaginary part (blue) of the prefactor 
term in the expression (|4.5p for A 3 (2, 2; a + it), with a G [3.5, 7.5], and t £ [0.1, 10]. 

phase values of the null contours there, again starting from zero on the real axis: 
-tt/2, -7t, 7t/2, 0, -7r/2, etc. 

We can use this information to understand the behaviour of the null contours 
shown in Fig. [3J As the null contours of Re(A 3 (2, 2m; s)) pass through a simple 
zero, their phase must change by n, from 7r/2 to —tt/2, or vice versa. However, we 
note that the relevant contours on the left in Fig. [3J which are almost symmetric 
to those on the right have opposite phase values. This means the null contours of 
Re(A3(2, 2m; s)) vary smoothly as they cross the critical line. For the null contours 
of Im(As(2, 2m; s)) the situation is different: the phase change of it forces them 
to avoid their almost symmetric counterpart on the left and "jump up" a contour 
as they pass through a zero on the critical line. These remarks are in accord with 
the derivative estimates at the critical line (see equation (|3 . 15[) and subsequent 
discussion). 

5. Some properties of zeros 

The properties of the angular sums given in the preceding sections enable simple 
proofs to be given of some properties concerning the location of zeros of angular 
lattice sums. These properties form an interesting counterpoint to those of the 
zeros of Epstein zeta functions (see Bogomolny and Lebouef (1994) and references 
[26-30] of that paper). The Epstein zeta functions are characterized as having an 
infinite number of zeros lying on the critical line, but with many zeros lying off 
that line and with almost all zeros lying on the critical line or in its immediate 
neighbourhood. We have already seen in the example of Section 2 that this third 
property of Epstein zeta function zeros is not shared by those of the trigonometric 
double sum C(2m, 1; s) for large m. 

Theorem 5.1. The trigonometric sums C(2,2m;s) and 5(2, 2m; s) have no zeros 
for s on the critical line in the asymptotic region t >> 4m 2 which are not zeros of 
both, and these will then be simultaneous zeros of C(0, 1; s) and (7(1, 4m; s). 
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Proof. We combine equations Q3.4p and (13. 5|) , to obtain 



C(2, 2m; s) + 5(2, 2m; ,s) = 



\/^2m(s) 



T(s) 



C(0,l;s) 



(5.1) 



As the term in square brackets is real on the critical line, from the functional 
equation (|3.ip . we see that 



arg[C(2, 2m; + S(2, 2m; ^ + it)] 



-4>2m;c(t) + 



(5.2) 



The second term on the right-hand side of (|5.2j) is chosen according to the sign of 
the term in square brackets in (|5.ip . We can also use (13. 9|) to deduce that 



arg[C(2, 2m; - + it) - 5(2, 2m; - + it)] 



(5.3) 



Note that the term on the right-hand side in (|5.3|) incorporates the phase of the 
real function on the left-hand side, and is not necessarily the same as the second 
term on the right-hand side of 



If we suppose that 5(2, 2r 



it) = 0, we require from (|5.2[) and (|5.3p that 



4>2m;c(t) = n7T f° r some integer n. We know from (|3 . 14[) that this cannot occur. 
Exactly the same argument applies if C(2, 2m; \ + it) — 0, so neither function can 
be separately zero on the critical line in the asymptotic region. If both are zero 
then, trivially, C(0, 1; s) = and C(l, 4m; s) = 0. □ 

We now consider the phases of 5(2, 2m; \ +it) and C(2, 2m; \ + it), particularly 
in the neighbourhood of zeros C(0, l;s) or C(l,4m;s). For brevity, we will adopt 
the following notations: 



\C(2, 2m; ~ + it)] = ]C\, argC(2, 2m; i + it) = 6 C , 



and 



|S(2,2m;i + a)| = arg5(2, 2m; 1 + it) 



it) = |5|(cos6 s 



Then on the critical line, we have for C and S that 
(7(2, 2m; - + it) = ]C] (cos 8 C + i sin 6 C ), 5(2, 2m; - 

and 

|C|sin9 c = |5|sin9 s . 

Using (|5.6p and (|5.7p . we can rewrite (|5.2p as 

2 sin C sin 9 S . . 

■ , n , n ; = -tan^ 2mjC (t). 
sm(B c + Q s ) 

With a prime denoting t derivatives, we obtain by differentiating ((57 

(2 cos 6 C sin e s )9' c + (2 sin 9 C cos Q s )0' s 



(5.4) 
(5.5) 



isine s ), 
(5.6) 



(5.7) 
(5.8) 



COS 2 </>2m,c(i) 



sin(6 s + C ) - tan (fe m , c (t) cos(8 ;j + 6 c )(e' s + 6^ 



(5.9) 
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Figure 7. Variation of the quantities <j>2,c{t) (red), c (i) (blue) and & s (t) (green) with 
t G [9.5,20] (left) and [20,30.5] (right), where those angles pertain to (7(2, 2; s) and 
5(2,2; s). 

We now consider the values of @ s , C and their t derivatives in the neighbour- 
hood of zeros of C(l,4m; \ + it) or C(0, 1; ^ + it), which are not zeros of both 
functions. Close to a zero of C(l, 4m; | + if), we require 



|c|c OS e c -|s|cose s ->o, e c ^e s , 



(5.10) 



where \C\ ^ 0, |5| at the zero by Theorem 5.1. Using the last of (|5~T0| in ([5 
and ()5.9p . we find at the zero of C(l, 4m; | + it) that 



, e' c + e' s = -2^ m , c (t). 



<d s = 6 C = -(t>2m,c(t) " 

Close to a zero of C(0, 1; \ + it), we require 
|C|cos6 c + 1^1 cos 8 S -> 0, \C\ sin 6 C + |5|sin0 s 0, \C\ -> 
From flUTH) , at the zero 

el = -e' 



(5.11) 



(5.12) 



(5.13) 



The relationships are exemplified in Fig. [71 where the variations with t of the 
functions cot(0 c ) and oot(0 s ) are shown for the case m = 1. Note that, in terms of 
these, the relationship (|5.8j) controlling the variation of A3 (2, 2m; s) on the critical 
line is 

COt(6 c ) + C0t(6 s ) = -2 COt(02m,e(t). (5.14) 

Then, from (j5TTTj) . the zeros of C(l, 4m; ^ + it) are evident where the curves giving 
cot(0 c ) and cot(0 s ) intersect, which they do on the line giving the variation of 
cot((j} 2m ,c(t)- The curves of cot(6 c ) and cot(6 s ) also go off to infinity at zeros 
of C(0, 1; \ + it), with their movement being in opposite directions in accord with 
(|5.14[) . In Fig.[JJwe see examples of there being none, one or two zeros (7(1, 4m; h+it) 
between successive zeros of (7(0, 1; \ + it). 

The next result relates to the independence of the functions A3 (2, 2m; s) on the 
critical line, again in the asymptotic region. 
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Theorem 5.2. If s is on the critical line in the asymptotic region t >> 4l 2 and 
ati, are non-zero reals, then 

aiC(l,Al; s) + 0C(0, 1; s) ^ 0, unless C(l, 41; s) = = (7(0, 1; s). (5.15) 

// s is on the critical line in the asymptotic region t >> 4max(( 2 ,m 2 ) and cti, a m 
are non-zero reals with I ^ m, then 

a t C(l, 41; s) + a m C(l, 4m; s) f unless C(l, 41; s) = = C(l, 4m; s). (5.16) 

Proof. We suppose 

a;C(l, 4/; s) + /?C(0, 1; s) = with (7(0, 1; s) ^ 0. (5.17) 

Then 

r.(\ ai- B \ \r,(\ ai- B \l 

(5.18) 



C(1,4Z; S ) 
O(0, l;s) 



C(l,4Z;a) 
0(0,1; a) 



for m an integer. 

We now use the functional equations (13. 1|) for C(0, 1; s) and C(l,4m; s), which 
on the critical line with 1 — s = s enable their arguments to be deduced: 

arg(7(0,l;s) = -[arg7r s + argT(s)] , argC(l,4m;s) = -[arg7r s + argT(s + 2m)}. 

(5-19) 

Using (|5.19p we arrive at the estimate for s on the critical line in the asymptotic 
region 

2m 2 1 

argC(0, l;s) - argC(l,4m; s) = mir + 0{ — ). (5.20) 

t t z 

Using ([5^0]) . (|5TT5|) requires 

2Z 2 1 

_ + O(-) - Ztt = m7r , (5.21) 

which is not possible. This proves (|5.15[) . 
Consider next (|5.16p . and suppose 

aiC(l, 41; s) + a m C(l, 4m; s) = with C(l, 4m; s) ^ 0. (5.22) 

Then this requires 



arg 



" C{l,4l;s) 
_C(l,4m;s) 

for an integer p. However, from (|5.20[) . 



Ptt, (5.23) 



2<l 2 - m 2 ) 1 

argC(l,4Z;s)-argC(l,4m;s) = (m-iW+J i + o( ). (5.24) 

Hence (|5.22j) requires 

2(P ; m2) + O(l) = (p + Z - m)7r, (5.25) 
which is not possible if / 7^ m. This proves (|5.16[) . □ 
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Corollary 5.3. The sums C{2, 2, s), 5(2, 2, s) and C(l, 4; s) have no zeros s on the 
critical line in the asymptotic region t >> 4 which are not zeros of both C(0, 1; s) 
and C(l, 4; s). 

Proof. We apply Theorem 5.2 and equation (I5.15|) . Now the sums C(2, 2, s), 5(2, 2, s) 
and C(l, 4; s) belong to the system of order 4 discussed in Appendix A, which is gen- 
erated by linear combinations of C(0, 1; s) and C(l, 4; s). All such combinations are 
of the form (|5.15[) , and are thus zero only if both C(0, 1; s) = and C(l, 4; s) = 0. □ 

We can also simply establish a corresponding result for the derivative A3 (2, 2m; s) 
for s on the critical line. 

Theorem 5.4. The derivative function A 3 (2, 2m; s) has no zeros for s = 1/2 + it 
in the asymptotic region t >> Am 2 which are not zeros of A3 (2, 2m; s) of multiple 
order. 

Proof. We recall the result O for s = 1/2 + it: 

A 3 (2, 2m; 1/2+it) = [1 -itan(0 2m . c (i))][|C(2, 2m; l/2 + i£)| 2 — |5(2, 2m; 1/2 + it)] 2 ], 

(5.26) 

where </>2m,c(i)) is a real-valued function. Taking the derivative of A3 (2, 2m; 1/2+it) 
with respect to 1/2 + it and separating into real and imaginary parts, we find 

A 3 (2, 2m; 1/2 + it) = - ( ^ 2m; 1/2 + it)\ 2 - |5(2, 2m; 1/2 + it)\ 2 ] 

+ tan0 2m , c (t)^[|C(2, 2m; 1/2 + zi) | 2 - |5(2, 2m; 1/2 + it)\ 2 ]^ 

-i-^[\C{2, 2m; 1/2 + it)\ 2 - \S{2, 2m; 1/2 + it)\ 2 ]. (5.27) 
For A3 (2, 2m; 1/2 + it) to be zero, we find from the imaginary part of (|5.27| that 

4[|C(2, 2m; 1/2 + it)\ 2 - 15(2, 2m; 1/2 + it)\ 2 ] = 0. (5.28) 
at 

We use (|5.28j) in the real part of (|5.27p . and, assuming tan^2m,c(^) is non-singular 
and <p' 2m Jt) 7^ 0, we complement (|5.28j) with 

[\C(2, 2m; 1/2 + it)\ 2 - |«S(2, 2m; 1/2 + it)\ 2 ] = 0. (5.29) 

Since both the assumptions we have just mentioned are true in the asymptotic 
region, we have proved the theorem. □ 

We provide plots of A 3 (2, 2; 1/2 + it) in Fig. [5] These show that the derivative 
is non-zero in the non-asymptotic region as well as in the asymptotic region. The 
numerical value of A 3 (2, 2; 1/2) is 0.918604. 

We next consider the properties of trajectories of constant phase which follow 
from the assumption that the Riemann Hypothesis holds for A3(2,2m;s). As can 
be seen from Fig. [HI lines of constant phase between zeros of A 3 (2, 2; s) do not in 
general cross the critical line, but asymptote towards it, and their configuration 
is arranged about a zero of the derivative <9|A 3 (2, 2; a + it)\/dt or equivalently of 
<9argA 3 (2,2;cr + ii)/dcr. 
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Figure 9. Contours of constant phase for As(2,2;s) in the region around its zero at 
s = 0.5 + il9.80599, with a second zero evident at s = 0.5 + 621.02204. The plot at right 
shows detail around the zero of the derivative <91og |Aa(2, 2; a + it)\/dt. 

Theorem 5.5. Given the Riemann hypothesis holds for A 3 (2,2m;s), then lines 
of constant phase coming from a = — oo can only cut the critical line at a zero of 
Aa(2, 2m;s) or of <9arg A^(2, 2; a + it) /da. Those lines of constant phase coming 
from a — oo cutting the critical line at a point which is not a zero of A3 (2, 2m; s) or 
of <9arg A3 (2, 2; a + it) /da must curve back and pass through a zero of A3 (2, 2m; s). 

Proof. The assumption of the Riemann hypothesis holding enables us to say that 
lines of constant phase coming from a = — 00 do not intersect before the critical 
line. Their phase monotonically decreases as t increases in a < 1/2, while the 
phase of A(2, 2m; 1/2 + it) monotonically increases with t. Thus, groups of constant 
phase lines coming in from a — —00 cannot cut the critical line, except at a zero 
of A(2, 2m; s). Isolated trajectories passing through points where <9arg A(2, 2; a + 
it) I da = are allowed. Such special trajectories separate lines of constant phase 
which curve up as a — > 1/2 and they move towards a zero of A(2,2m;s), from 
trajectories which curve down as a — > 1/2. 

Now consider sets of lines of constant phase approaching a =1/2 from the right. 
Their phase increases with t, but if they were able to cross the line a = 1/2 and 
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progress towards a — — oo it would have to decrease with t. Thus they must turn and 
run alongside a = 1/2, with lines above the special trajectory curving up towards 
a zero of A(2, 2m; s), and those below curving down towards a zero. Alternatively, 
they can cut the critical line at a point which is not a zero of A 3 (2, 2m; s) or of 
9argA(2, 2; a + it) /da. Such trajectories supply the lines of constant argument 
required for generic points on the critical line, and must return back to cut the 
critical line at a zero of A(2,2m;s). The region where they cross into a < 1/2 is 
bounded on the left by a line of constant phase connecting a zero of A3 (2, 2m; s) 
with an adjacent zero of d arg A(2, 2; a + it)/ da. □ 

Note that the region where zeros of 9 arg A (2, 2; a + it) /da exist is confined 
to the neighbourhood of the critical line, since this partial derivative being zero 
corresponds to a horizontal segment on a line of constant phase (see, for example, 
Fig. [3]). However, we may use the line of constant phase passing through the point 
on the critical line where d arg A (2, 2; a + it) /da = as the separator between lines 
of constant phase going to the zero above this line from those going to the zero 
below (given the Riemann hypothesis is assumed to hold). 

Theorem 5.6. Given the Riemann hypothesis holds for A3(2,2m;s), then there 
exists one and only one zero of d arg A3 (2, 2; a + it) /da on the critical line between 
two successive zeros of A3 (2, 2m; s). 

Proof. We consider the analytic function 



The real part of this function goes to —00 at any zero of A3 (2, 2m; s), and 
increases away from these logarithmic singularities. It must have at least one turning 
point between successive zeros. By the Cauchy-Riemann equations, such a turning 
point is a zero of d arg Aa(2, 2; a + it)/ da. 

Next, suppose there two or more zeros of 9 arg A(2, 2; a + it) /da between suc- 
cessive zeros so = 1/2 + Uq and si = 1/2 + it\. Denote the upper two of such 
derivative zeros by s* = 1/2 + it* and s** = 1/2 + it**. As we have seen, each of 
these has constant phase lines coming from a = — 00 and passing through it, around 
which constant phase trajectories reverse their course. Those above s* curve up to 
si as they approach the critical line, while those below it curve down. They cannot 
cross the constant phase line passing through s**, nor can they cross the critical 
line. They must then head back to a — —00, where they will breach the mono- 
tonic nature of the variation of arg A3 (2, 2m; s) with t. Thus, this situation cannot 
arise. □ 

We conclude this section with an investigation of the structure of lines of con- 
stant phase which cut the critical line and turn back to a = 00 thereafter. We start 
with the expansion of arg A3 (2, 2; a + it) around a point s* = a* + it* = 1/2 + it* 
in the asymptotic region where d arg A3 (2, 2; a + it) /da — 0, which is of the form: 



logA 3 (2,2m;s) = log | A 3 (2, 2m; s)| + i arg A 3 (2, 2m; s). 



(5.30) 



arg A 3 (2, 2; o- + «t) 



arg A a (2, 2; 1/2 + it.) + 




(a-1/2) 2 



2m 2 



{t-t*f 



2mr T d 2 

+ ^2log|A 3 (2,2;a + i t)| 

(a - a*)(t - t*) + . . . . 



(5.31) 
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0.4995 0.4996 0.4997 0.4998 0.4999 0.5000 0.5001 0.4992 0.4994 0.4996 0.4998 0.5000 



Figure 10. Contours of constant phase for As(2,2;s) in the region around its zero at 
s — 0.5 + H9. 80599. The thick line is the contour on which <9arg As(2, 2;a + it)/dt is zero. 
(a)The fine lines correspond to the phases 3.04,3.04076, 3.0410, 3.0412, 3.0414, 3.0416, 
3.0418, 3.0419, 3.04198, 3.042, 3.0425, 3.043, 3.044, 3.045, 3.046, 3.047 and 3.048. (b) The 
thick dashed line is the contour on which 9 log |Aa(2, 2; a + it)\/dt is zero. The contours 
are -0.1090, -0.1080, -0.1072, -0.1068, -0.1064, -0.1060, -0.1056, -0.1055, -0.1054, -0.1053, 
-0.1052, -0.1048, -0.104, -0.103, -0.102, -0.1015 and -0.10. 



Here we have employed the asymptotic estimates for arg Aa(2,2; 1/2 + it) based 
on (|3.14| . The corresponding trajectories of constant phase may be shown to be 
rectangular hyperbolae, with their centre at 



(7 - 1/2 



t 



l[^log\A 3 (2,2;o- + it)\] s=st 
1 



4 L('/ 



4m 4 



2tl 
{2m 4 



d2 log |A 3 (2, 2; a + ii)!]^, V H 



dt 2 



log 1^(2,2; (7 + it)| 



(5.32) 



The second derivative factor in (|5.32p is always negative, and tends to be much 
larger than the terms involving powers of 1/i*. Thus, we see that the centre of the 
hyperbolic trajectories of constant phase will always lie to the left of the critical 
line, with its ordinate very close to i*. This displacement of the centre into a < 1/2 
creates the region in which phase lines can cut through a = 1/2 and return to a > 
1/2 via passage through a zero of As(2, 2; s). Note that the centre corresponds to a 
point at which two lines of constant phase intersect; thus, it must have derivatives 
of phase along two independent lines which are zero. It therefore is a point at 
which both darg A 3 (2, 2;a + it)/dt and <9arg A 3 (2, 2; a + it) /da are zero. From the 
Cauchy-Riemann equations, it also is a point at which 91og |A3(2, 2; a + it)\/dt and 
<91og |Aa(2, 2; a + it)\/da are zero- i.e., it is a point of extremum for both amplitude 
and phase. 

These hyperbolic centre points are locations at which A3 (2, 2; s)' = 0. The fact 
that they lie to the left of di(s) = 1/2 is is similar to the property proved by Speiser 
(1934) for C(s). 
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The curves of constant phase of A3(2, 2; s) given in Fig. ITD1 illustrate some more 
of their general characteristics in the neighbourhood of zeros. For Fig. [10] (a), at 
the top of the figure, we see the centre of the hyperbolic phase curves; through this 
centre passes a curve on which darg A 3 (2, 2; a + it)/dt = 0. This curve connects 
the hyperbolic centre to the zero of A 3 (2,2;s) below it, and also continues to the 
zero above; at each zero it is tangent to the t axis. The curve marks points where 
the constant phase lines have vertical slope. The lowest line of constant phase on 
the right has vertical slope when it arrives at the zero of Aa(2, 2; s); it corresponds 
to 3.04076, which is the value given by (|3.8|) for the phase on the critical line 
just above the zero. A second important line has the phase 3.04198, which is the 
phase corresponding to the centre of the hyperbola. This line again connects the 
hyperbolic centre to the zero. All lines whose phase lies between these values come 
in from the right, cross the zero line of darg A3 (2, 2; a + it)/dt, and curve back to 
pass through the zero of A 3 (2,2;s). Where they lie to the left of the zero line of 
darg A 3 (2, 2; a + it)/dt, their phase increases as t decreases; where they lie to the 
right, it increases as t increases. Curves coming in from the left all pass through 
the zero without crossing the zero derivative line; their phase always increases as t 
decreases. 

For Fig. [10] (b), we show the phase contours below the zero, in the region down 
to the next zero. The phase at the centre of the hyperbolae is -0.1055, while the 
phase just above the bottom zero and just below the upper zero are respectively 
-0.107604 and -0.101134. Note the dashed line passing through the centre, along 
which 9 log |A3(2, 2; a + it)\/dt = 0. From the Cauchy-Riemann equations, this 
is perpendicular to the solid line defined by darg A3 (2, 2; a + it)/dt = (a fact 
disguised by the different scales on the horizontal and vertical axes). 



6. Distributions of zeros 

We return to the left-hand side of expression 1|4.6|) . in which we replace s by a + it, 
and expand assuming \t\ >> a, with a large enough to ensure accuracy of (|4.6p . 
The result is 

2ilog(f) - 2i(log7r + 1) + 7r(er - -) + j(l - 2a). (6.1) 

As for each increment of tt of this expression we get one null line of the real part 
of A3 and one of the imaginary part, and these intersect at a = 1/2 to give one 
zero there (assuming the Riemann hypothesis holds for A3 (2, 2m; s)), we can divide 
(|6.ip by tt, and regard the result as a distribution function for zeros of A3: 

iV A ,(*, t) = - log(i) - -(1 + logTr) + a - ~ + ^(1 - 2a). (6.2) 

Now, from Titchmarsh and Heath-Brown (1987), the distribution function for 
the zeros of the Riemann zeta function on the critical line is 

N d\,t) = £ Iog(*) - ^(1 + log(27r)) + O(logi). (6.3) 

We complement this with the numerical estimate from McPhedran et al (2007) for 
the distribution function of the zeros of L_4(s): 

N-A,t) = ^log(t) - ±(1 + log(7r/2)) + O(logt). (6.4) 
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1 2 3 4 5 1 2 3 4 



Figure 11. The distributions of the differences between successive zeros of C(0, 1; s)(left) 
and C(l, 4; s) (right) on s = 1/2 + it, for i < 300. 

Adding (|6.3[) and (|6.4[) we obtain the distribution function for the zeros of C(0, 1; s) 
(see (1531) ): 

JVco,i(i*) = - log(t) - -(1 + log(Tr)) + O(logt). (6.5) 

Z 7T 7T 

When we compare this with (|6.2|) , and use the equation 

N A3 (\,t) = N C1A (^ t) + N C0 ,i(~,t), (6.6) 

it suggests the hypothesis that the distribution function of zeros of C(l, 4; s) is the 
same as that of (|6.5[) , to the number of terms quoted: 

JVci >4 (~,t) = Afco.id. t) = - log(<) - -(1 + log(Tr)) + O(logi). (6.7) 

l Z 7T 7T 

Strong numerical evidence supporting this is given in Table [TJ which also shows 
zero counts for C(l, 8; s) and C(l, 12; s). Note that the numbers of zeros found for 
C(l, 4; s), C(l, 8; s) and C(l, 12; s) are virtually the same. This rules out any variation 
with increasing order similar to that of Dirichlct L functions, where increasing order 
results in significant increases in density of zeros (compare (|6.4|) and (|6.3|) . or the 
second and third columns of Table 1). 

Comparing the data of Table 1 with the discussion in Bogomolny and Lebouef 
(1994), we can see that the split up of N(t) into averaged parts given by expres- 
sions like (j6.3H6.7p and oscillating parts applies to C(0, 1; s) and to the C(l, 4m; s). 
However, it would be value to extend the numerical investigations of Table 1 to 
much higher values of t, to render the characterization of the oscillating term more 
accurate. Such an extension may require the development of an alternative algo- 
rithm to that based on (|2.3p . which will probably become unwieldy for values of t 
of order 10 4 - 10 5 . 

In Fig. [Tl] we compare the distributions of the differences between zeros on 
the critical line for C(0, 1; s)(left) and C(l, 4; s) (right). The distributions are quite 
different, even with this modest data set. Bogomolny and Lebouef (1994) have 
studied the case of C(0, l;s) using 10,000 zeros after t = 10 5 , and contrast the 
distribution for £(s)L_4(s) with that for each function separately. The separate 
factors in fact have distributions like that that for C(l, 4; s). The function compared 
with the histogram in the right of Fig. [Tl] corresponds to the Wigner surmise, 
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which (Bogomolny and Lebouef, 1994, Dietz and Zyczkowski, 1991) for the unitary 
ensemble takes the form 



Here the separation between zeros has been rescaled to have a mean of one. Bo- 
gomolny and Lebouef (1994) comment that the left distribution is that of an un- 
corrected superposition of two unitary ensemble sets. Note as one indicator of this 
that there is not the same pronounced tendency for the probability to go to zero 
with separation on the left as on the right, where the distribution clearly comes 
from a single ensemble. Recall that Table 1 shows that the frequency distribution 
for zeros is the same for C(0, 1; s) and C(l,4; s). This makes the strong difference 
in the distributions of the gaps all the more interesting. 

The research of R.McP. on this project was supported by the Australian Research Council 
Discovery Grants Scheme. 



Abramowitz, M. & Stegun, I. A. 1972 Handbook of Mathematical Functions with Formulas, 

Graphs and Mathematical Tables. New York: Dover. 
Bogomolny, E. & Lebouef, P. 1994 Statistical properties of the zeros of zeta functions- 

beyond the Riemann case. Nonlineanty 7 1155-1167. 
Dietz, B. & Zyczkowski, K. 1991 Level-spacing distributions beyond the Wigner surmise. 

Z. Phys. B- Condensed Matter 84 157-158. 
Gonek, S., 2007 Finite Euler Products and the Riemann Hypothesis arXiv:0704.3448 
Hardy, G. H. 1920 On some definite integral considered by Mellin. Mess. Math. 49, 86-91. 
Kober, H. 1936 Transformation formula of certain Bessel series, with reference to zeta 

functions Math. Zeitschrift 39, 609-624 
Lorenz, L. 1871 Bidrag til talienes theori, Tidsskrift for Math. 1 , 97-114. 
McPhedran, R. C, Smith, G. H, Nicorovici, N. A. & Botten, L. C. 2004 Distributive and 

analytic properties of lattice sums. J. Math. Phys. 45, 2560-2578. 
McPhedran, R. C, Botten, L. C. & Nicorovici, N. A. 2007 Null trajectories for the sym- 
metrized Hurwitz Zeta Function Proc. Roy. Soc. A 463 303-319 
McPhedran, R.C, Zucker, I.J., Botten, L.C. & Nicorovici, N.A.. 2008 On the Riemann 

Property of Angular Lattice Sums and the One-Dimensional Limit of Two- Dimensional 

Lattice Sums. Proc. Roy. Soc. A, 464, 3327-3352. 
Speiser, A., 1934 Geometrisches zur Riemannschen Zetafunktion. Math. Ann. 110, 514- 



Titchmarsh, E. C. & Heath-Brown, D. R. 1987 The theory of the Riemann zeta function, 

Oxford: Science Publications. 
Zucker, I.J. & Robertson, M.M.. 1976 Some properties of Dirichlet L series. J. Phys. A. 

Math. Gen. 9 1207-1214 



Article submitted to Royal Society 




(6.8) 



References 



521. 



New Properties of Angular Lattice Sums 



21 



Table 1. Numbers of zeros 0/ C(l/2 + it), L_ 4 (l/2 + it), C(l, 4; 1/2 + \t),C{l, 8; 1/2 + it) 
and C(l, 12; 1/2 + it) in successive intervals of t. 







t 


n c 


n_4 


ncu 


n c + 
n-4 + ncn 


nci8 


WC112 


0-10 





1 


2 


3 


2 


3 


10-20 


1 


4 


5 


10 


5 


5 


20-30 


2 


5 


6 


13 


7 


7 


30-40 


3 


4 


8 


15 


8 


8 


40-50 


4 


6 


8 


18 


8 


8 


50-60 


3 


5 


9 


17 


9 


9 


60-70 


4 


6 


9 


19 


10 


10 


70-80 


4 


6 


11 


21 


10 


10 


80-90 


4 


7 


11 


22 


11 


10 


90-100 


4 


6 


10 


20 


10 


12 


0-100 


29 


50 


79 


158 


80 


82 


(P^l).(t6^.(l67)l 


28 


50 


78 


156 






100-110 


4 


7 


11 


22 


11 


10 


110-120 


5 


7 


12 


24 


12 


12 


120-130 


5 


7 


12 


24 


12 


12 


130-140 


5 


7 


12 


24 


12 


12 


140-150 


4 


7 


12 


23 


12 


11 


150-160 


6 


7 


12 


25 


12 


13 


160-170 


6 


7 


13 


26 


13 


12 


170-180 


6 


8 


13 


27 


13 


13 


180-190 


5 


8 


13 


26 


12 


14 


190-200 


5 


7 


13 


25 


14 


13 


0-200 


80 


122 


202 


404 


203 


204 


(|Ql).(p^)l.(f6T7l) 


78 


122 


200 


400 






200-210 


6 


8 


13 


27 


12 


13 


210-220 


5 


8 


14 


27 


15 


14 


220-230 


6 


8 


13 


27 


13 


13 


230-240 


6 


8 


14 


28 


14 


14 


240-250 


6 


8 


14 


28 


14 


13 


250-260 


6 


8 


13 


27 


14 


15 


260-270 


6 


8 


15 


29 


14 


13 


270-280 


6 


8 


14 


28 


14 


15 


280-290 


6 


8 


15 


29 


15 


14 


290-300 


6 


9 


14 


29 


13 


15 


0-300 


137 


203 


341 


681 


341 


343 


(p^l). (^4)1.(1671) 


137 


203 


340 


680 







Appendix A. Supplementary Notes on Systems of Angular 

Sums 

We give here results linking trigonometric lattice sums of order up to ten, which 
show that they may be generated from three independent sums, C(0, 1; s), C(l, 4; s), 
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and C(l, 8; s). We also give expressions in terms of sums of Macdonald functions of 
the Kober-type from which these independent sums may be calculated. 
A basic result we use relies on the symmetry of the square lattice: 



C(l, 4m - 2; s) = T > cos ^™ Wvuv* = q = c(2j 2m _ 1}s) _ 5(2j 2m _ 1; s)) 

PUP2 

from which we obtain 



C(2,2m- l;s) =5(2, 2m- l;s) = -C(0,l;s), (A2) 

for m any positive integer. We can also expand cos(4m— 2)0 PliP2 = T4 m _2(cos # PliP2 ) 
in terms of powers of cos Pl jP2 , using the expressions for the Chebyshev polynomials 
in Table 22.3 of Abramowitz & Stegun (1972). This enables us to express C(4m — 
2, 1; s) in terms of C(in, 1; s) with An < Am — 2, and n being a positive integer. In 
this way, we can inductively arrive at the results below. 

(a) Order 

The only sum of this type is C(0, 1; s), given by equation (|2.2|) . It is given by a 
modified form of (|2.3j) . since there are contributions from both axes p% = 0, p2 = 
rather than just p2 = 0: 

C (0,l; S ) = 2C(2 S ) + 2 ^;- 1 / 2 ) C (2,-l) 

r(s) 

„ „ oo oo / \ s— 1/2 



E E (jf) ^-iM 2 ^)' (A3) 



r(s) ^ ^ \pi 

pi = l P2=l 

(6) Order 2 



The single sum of order 2 is 



£ W^ = C(2 ' 1;s) = 5 C( °' 1;s) 

(Pl,P2) 

i^+i) + i>Ti) ? ?U J m* k *iM**p#>)- 



(A 4) 



(c) Order 4 

We generate this system from C(4, 1; s), obtained from (|2.3p . In terms of this, 
C(l, 4; s) = 8C(4, 1; s) - 3C(0, 1; s), (A 5) 

and 

C(2,2;s) = 4C(4, l;s) -C(0, l;s), 5(2, 2; s) = -4C(4, 1; s) + 2C(0, 1; s). (A 6) 
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(d) Order 6 
The result of expanding (|A ip for m = 2 is 



while from (IA 2 



Also, 



C(6,l;«) = |c(4,l;«)-ic(0,l;«), (A 7) 



C(2,3;«)=5(2,3;«) = ic(0 ) l;«). (A 8) 



^ (Pi+Pl) s+3 ^ (Pl+P 2 2 ) s+3 2 ^ (pf+p|) s + 2 ' 1 j 

(Pi.pa) (Pi,P2) (pi,P2) 



so 



\ - / Cos 4 (^Pi,P2)sin 2 (6 | p 1 ,p 2 ) _ 1 , > . 
(P1.P2) 

Note that all sums mentioned so far can be generated from just two, say C(0, 1; s) 
and C(l,4;a). C(6, l;s) is given by (|2~3| . 



(e) Order S 

The new sum we use here is C(8, 1; s), obtained from (|2.3p . In terms of this, 

C(l, 8; s) = 128C(8, 1; s) - 224C(4, 1; s) + 49C(0, 1; s), (A 11) 

C(2, 4; s) = 64C(8, 1; s) - 112C(4, 1; s) + 25C(0, 1; s), (A 12) 

and 

5(2, 4; s) = -64C(8, 1; s) + 112C(4, 1; s) - 24C(0, 1; s). (A 13) 

Two other sums in this system are 

E \/3 Y+i = "C(8, 1; a) + §C(4, 1; s) - \c(0, 1; S ), (A 14) 

and 

E ' ? 2 I lP L +4 - c ( s > !; s ) - 2C ( 4 > !; s ) + o c (°' !; «)■ ( A15 ) 

( Pl ,p 2 ) U'i+w 

Using the result (|A lip , we have calculated curves showing the modulus of 
C(l,8;s) as a function of s = 1/2 + it on the critical line. These are given in 
Figs. fl2lfT3l and the distributions of zeros they show are given in Table [TJ 
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Figure 12. The modulus ol C(l, 4; s) (red) and C(l, 8; s) (blue) as a function of 
s = 1/2 + it, for t SE [0, 20] (left) and t € [20, 40] (right). 




Figure 13. The modulus of C(l, 4; s) (red) and C(l, 8; s) (blue) as a function of 
s = 1/2 + it, for t € [40, 60] (left) and t € [60, 80] (right). 



(/) Recurrence Relations 

We can generalize the above procedure by establishing recurrence relations for 
the trigonometric sums. We consider 



2n 
Pi 



(P1+P2) 



2^s+n 



= C(2n,l;«)= ^ 



1 



IPi + Pa) 



(P1.P2) ^ (P1,P2)^(0,0) 

Expanding using the Binomial Theorem, we obtain 

n 

C(2n, 1; s) = ^ n Q C(2l, 1; s). 



1 



Pi 



Pi +vl 



1=0 



If n is even, (|A 17[) gives an identity: 



2n-l 



£ a «c I (-i) , c(2Z,i;») = o ) 



(A 16) 
(A 17) 

(A 18) 



2=0 



while for n odd we obtain an expression for C(4n — 2, 1; s) in terms of lower order 
sums: 



2n-2 



1 - 

C(4n- 2, !;«) = - £ ^C, (-1)^(21,1; a) 



(A 19) 
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5 10 



Figure 14. Null contours ol the real part (red) and imaginary part (blue) of 
A 3 (2,4;cr +it), with a G [-4.5, 11.5], and t SE [0.1,20]. 



It may be checked that the two relations (|A 18[) and (|A 19[) give equivalent results 
in the cases given above. 

In Fig.[l4]we give null contours of the real part and imaginary part of As(2, 4; s), 
for comparison with those of Aa(2,2;s) in Fig. [3] It will be noted that there are 
now three lines giving the contours on which Im A3 (2, 4; s) = 0, and these are not 
circles (as was the single exemplar in Fig. [3]). There is a corresponding increase 
in the number of null contours starting and finishing on the real axis, and in the 
value of a at which the last one reaches the real axis. There are two examples 
in Fig. [14] of zeros on the critical line of the real part, but not the imaginary part 
(tan (02m,c(l/2 + it) = 00), and one of a zero of the imaginary part but not the real 
part (cot (</>2m,c(l/2 + it) = 00). The null contours starting at a = 00 and ending 
at a = — 00 settle down to an asymptotic behaviour for somewhat larger values of 
t than in Fig. 02 but are otherwise similar to the lower order null contours. 
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(g) Order 10 

The new sum we use from (|2 .3[) is C(10, 1; s), which, from (|A 19[) . is 

C(10, 1; s) = ~C(0, 1; a) - |c(4, 1; s) + lc% 1; s). (A 20) 



Other sums in the system are: 



V- / p\p\ 



r ^ (P?+Pl) s+5 

(P1,P2) 



C(8,l;s)-C(10,l;s), (A 21) 



and 



6 4 1 4 4 

, PlPl 1 / P1P2 



(PI1P2) (Pi,Pa) 



(A 22) 



Appendix B. Derivation of the Functional Equation for 

C(l,4m; s) 

We start with the equation (30) from McPhedran et al (2004). This equation uses 
the Poisson summation formula to derive a connection between a sum over a direct 
lattice of points R p = R pliP2 = d(p\,p2) with polar coordinates (R p ,4> p ), and a 
corresponding sum over the reciprocal lattice, where the lattice points are labelled 
K/j = K/ ll) ^ 2 . Considering the case of a square lattice with period d, the reciprocal 
lattice points are K/j = (27r/<i)(/ii, fe)- The sum in the direct lattice incorporates 
a phase term of the Bloch type, with wave vector ko, and the result of the Poisson 
formula is 

«2 S I,-,/ * n .e^ ^e^p 

(pi,P2) P 

27ri ; J . / e i/9 ° v-^ , e Mh \ 

= ^ r( 2 +1 - S) U^ S 'gpi ■ ( B1 ) 

Here the sum over the reciprocal lattice in fact runs over a set of displaced vectors: 

Q h = k + K h = (H^-l + fe 0:c , ^ + fcoy) = (Qh,0h), (B2) 



where in (|B 2ft the second and third expressions are in rectangular and polar coor- 
dinates. Note that in (|B 1|) and (|B 2|l . pi,P2 and h±,h,2 run over all integer values, 
and #0 gives the direction of ko. 

We express the relation (|B II) in non-dimcnsionalized form, taking out a factor 
d 2s on the left-hand side, and a factor (27r/ci) 2 ~ 2s on the right-hand side. We put 
ko = (2w/d)no, and obtain 

2 } ^ {pI+pD s V 2 

(Pl,P2) 
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If we now let kq — > 0, set I = Am and take Re(s) > 1, we obtain the desired result: 
QSLtf^,,^., - r ' 2 "^-» C(l,4m, 1 - .). (B4) 
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